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Abstract 

A scheme to reduce translational noninvariant quasi-one-dimensional wave guides 
into singly or multiply connected one-dimensional (ID) lines is proposed. It is meant 
to simplify the analysis of wave guides, with the low-energy properties of the guides 
preserved. Guides comprising uniform-cross-sectional sections and discontinuities 
such as bends and branching junctions are considered. The uniform sections are 
treated as ID lines, and the discontinuities are described by equations sets connect- 
ing the wave functions on the lines. The procedures to derive the equations and to 
solve reduced systems are illustrated by examples, and the scheme is found to apply 
when the discontinuities are distant and the energy is low. When the scheme applies, 
it may substantially simplify the analysis of a wave guide, and hence the scheme 
may find uses in the study of related problems, such as quantum wire networks. 

Key words: Wave guide, discontinuity, quantum wire, one-dimensional, 
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1 Introduction 

When a wave propagates with a low energy in a quasi-one-dimensional (QID) 
wave guide with an uniform cross section, the guide is effectively one-dimensional 
(ID) since only the first transverse mode plays a role. However, when a guide 
comprises discontinuities such as bends and branching junctions, higher trans- 
verse modes come into play at the vicinity of the discontinuities. Therefore, in 
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principle, a full-mode or full- wave analysis is required, which means a substan- 
tial increase in the amount of calculation and the physics is often obscured. 
Knowing that the guide is essentially ID in the uniform sections, such a full- 
wave analysis for the entire system is actually unnecessary. It is the purpose 
of this paper to present a scheme to eliminate those redundancies and simplify 
the calculations, yet have those low-energy properties faithfully preserved. 



Historically, there has been a number of schemes related to this purpose. The 
earliest one was due to Kuhn and dates back to 1949 [1]. Then it was made 
more well-known by Griffith in 1953 and henceforth it has been called the 
Griffith boundary condition [3|Hl5|6|7f8|9|10|lllll2lll3|14|15|16] . The scheme 
contains a set of equations relating the wave functions and their first deriva- 
tives on the lines connected to a junction [T|2|3|lf5|6|7|l8f9l[10|ll|12|13|14lll5lll6] 
The equations are simple equations that satisfy the unitarity condition at the 
junction. The scheme has been intuitively stated and has no undetermined 
parameters. However, recently it has been pointed out by the author et al. 
[T7] that it is not clear what kind of realistic QID guides the scheme de- 
scribes. Later than Kuhn, there was another scheme by Shapiro in 1983 [18]. 
The scheme starts by an unitary matrix relating the amplitudes of the inward 
and outward waves [TSP^EOj . Though this scheme contains free parameters, 
it also has not been mentioned how the parameters are related to the guides 
in realistic spaces. 



To the present, approaches to the related problems are belong to either one of 
the mentioned two categories [T|l^|5Pf5|5|7f51IMTU|lllll2 l!13!'14f 
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and a common feature of most of the schemes 
is that the relation between the reduced systems and the original systems is not 
addressed, and this hinders the application to the study of realistic systems. 



In 2006, the author et al. proposed a scheme [T7], which has been the first at- 
tempt to relate reduced and unreduced systems. The scheme resembles that by 
Kuhn |1], but has an extra phenomenological term with a tuning parameter. 
The parameter is to be fixed by a comparison between results from the reduced 
and unreduced systems. The scheme was shown to be a substantial improve- 
ment, and many low-energy transport properties were shown to be captured. 
But still, it is pointed out later in this paper that this simple phenomenology 
can be inadequate and a more general scheme is needed. 



Section [2] illustrates the derivation of the connecting equations to be used in 
reduced systems, for two typical component structures — the L-bend and the 
T-junction. Then in Sec. |3l these two structures are assembled into more 
elaborated structures, and the wave propagation in the structures are studied 
in and compared between, the reduced and unreduced systems. In Sec. IH a 
few concluding remarks are given. 
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2 Formulation 

The scheme to be proposed can be summarized as in the following. For a 
particular discontinuity, the scattering matrix (^-matrix) is evaluated using a 
full- wave treatment, and then the matrix is truncated leaving only those ele- 
ments relating the first transverse modes in the branch guides. The truncated 
S-matrix is then used to connect the ID wave functions in reduced systems. 

The two-dimensional (2D) time-independent Schrodinger equation (TISE) — [/i^/(2mo)][9|^'-|- 
dy"^] + V{x, y)"^ = E'if, with V{x, y) = and hard wall boundaries is consid- 
ered. For a wave with an energy E in a. guide having a width W, the wave 
function can be written as a sum of direct products of transverse modes and 
longitudinal waves. Labeling the guide by t] and defining a coordinate system 
(x^, yn) in the guide, where is in the longitudinal direction and is in the 
transverse direction (0 < < W), the wave function can be written as 

%{E;x,,y,)=Y. ^-sm [^-^ 

m=l I ' ' \ r r / 

where k^"^'^ = \j2mQE/h^ — (mTc/Wy is the longitudinal wave number for the 
m-th transverse mode, which can be propagating or evanescent, and iV^ is a 
large enough integer. In a ID space, the TISE becomes -[h^ / {2mo)]d'^iJ + 
V{x)il! = Eip. For a line labeled by i], and with a coordinate defined on it, 
the wave function for V{x) = is 

x„ y,) = A.e'''^- + B^e-'^^\ (2) 

where k is the longitudinal wave number given by = y/2mQE /h. When 
results from 2D guides and ID lines are compared, A^^^ (B^^^) is compared 
with (B,^), and k'^^^ is compared with k. 

Two discontinuities in 2D wave guides are considered, the L-bend and T- 
junction as shown in Fig. [TJ^a) and [11(b) respectively. The discontinuities are 
divided into regions which are labeled by 1, 2, 3, and D in our discussion. 
Coordinates are also defined in the branch guides, and wave functions in the 
branch guides are in the form given by Eq. ([T]). The S'-matrices are to be 
presented in terms of these coordinates. The symbols used to denote the L- 
bend and T-junction in reduced systems are shown in Fig. [^c) and [11(d) 
respectively. 

We only sketch the evaluation of the S'-matrices here, since the techniques are 
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well-established and are detailed in the literatures [2T|22] . The wave function 
in region D [see Figs. [U^a) and[T]^b)] is also expanded in terms of undetermined 
amplitudes, and they are connected to wave functions in other regions by the 
conditions of continuities of wave functions and normal derivatives of wave 
functions at the boundaries between the regions. Otherwise, one may also solve 
the TISE in a discretized space, where the TISE is a set of finite-difference 
(FD) equations. 

For the T-junction shown in Fig. [I](b), we can get an equation for the ampli- 
tudes of the waves in the branch guides (regions 1, 2, and 3) such as 
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where Aj'") (5^™)) is an amplitude for an inward (outward) wave in Eq. ([1]). 
Truncating the matrices in Eq. ([3]) and retaining only terms related to the 
first transverse modes in the branch guides, we get 
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where Ajj and i?^, r] = 1, 2, and 3 are amplitudes in Eq. ([2]), and 



(5) 



where the coordinates are chosen as toward the junction, and xi = X2 = = 
at the junction as shown in Fig. [T](d). For the scheme to work, it is necessary 
that the energy is such that only the first transverse mode is propagating, and 
the exponential tails of the evanescent waves of higher modes emanated from 
the discontinuities are shorter than the distances between the discontinuities. 

We can also write Eq. (jl]) in terms of iprj- Using = [ipj^ + dipri/d{ikXn)]/2 
and = [ipn — dipri/d{ikxrf)]/2 at = 0, the connecting equation can be 
rewritten as 
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where the wave functions and their derivatives are evaluated at the discon- 
tinuity, and the directions of the coordinates are defined to be toward the 
discontinuity. Note that in Eq. ([6]), it is not necessary that the origins of the 
coordinates be located at the discontinuity. 

Likewise, the connecting equation for the L-bend shown in Fig. [T](a) can be 
written as 



1 + Si 
ik 



dtjji 

dxi 
di/j2 

- dxo 



^1 

^2 



(7) 



where the directions of the coordinates are defined to be toward the disconti- 
nuity as shown in Fig. [H^c). 

Note that, the S-matrices are symmetric (^l* = 5'l and St^ = ) and 
unitary (Sl^^l = 1 and St^St = 1). The unitarity implies I]r?(l^?7p ~ 
l-Bj^n = 0, and since = [ip^ + d^prJ/d{ikxr,)]/[2e'^''^'^] and 5^ = [i/j^j — 
dil!ri/d{ikxr^)]/[2e~^^^'^], the unitarity can be rephrased as a more intuitive 
equality J2r] ^J^dipn/diikx^) = 0, or there is no net inflow of probability current 
to the discontinuity. 

The numerical results for the magnitudes and arguments of the elements of 
are plotted in Figs.[2](a) and[2](b) respectively, and the results for St are plotted 
in Figs. MyC) and El^d) respectively, versus a dimensionless longitudinal wave 
number k defined hy n = k^^^W/ii. For the reference of the readers, the cutoff 
of the second transverse mode is at k = a/3 ~ 1.73. The results are obtained 
using discretized spaces with 20 sites across a width of W, and they are found 
to be in congruence with results from 10 sites within windows of A/t ~ 0.02 on 
the horizontal axes, and windows of A\S\ ^ 0.02 and A[7r^^Arg(S')] ~ 0.005 
on the vertical axes, which implies that the continuous space limit has been 
approached. 

Since the L-bend and T-junction appear quite often in practical problems, it 
may be convenient to have their S'-matrices in analytic forms. Within a finite 
range of k, it is possible to approximate a 5'-matrix by analytic functions. For 
the S'l at K < 1, we may approximate the magnitudes and arguments of the 
elements by 



5 



IV Lyiii 1 + (2.796fi;- 1.498fi:2)2' ^ ' 

Arg (5l)ii ~ (1 - 0.586k + 1.496^2 - 0.541/t3)7r, (9) 

|(5l)i2| = V1-|(^l)iiP, (10) 

Arg(5L)i2 = Arg(5L)ii-^, (11) 

(5l)2i = (5l)i2, and (12) 

(5l)22 = (5l)ii. (13) 



Within K < 1, the approximation for K^l)!!] in Eq. ([8]) and the approximation 
for Arg (S'l)ii in Eq. ([9]) approximate the numerical results in Figs. [2]J^a) and 
[2]J^b) respectively up to ^\S\ < 0.01 and A[7r~^Arg(S')] < 0.01 accuracy. The 
exact equalities Eqs. ffTUl) - ffT^ are due to the symmetry and unitarity of the 
S'-matrix and the exchange of leads 1 and 2. 

For the St at k < 1, it is found that the numerical results for the elements 
can be approximated by 



+ (1.734.- 0.808.T' ^'^^ 

Arg (5t)ii ^ (1 - 0.051k + 0.559.2 ^ 0.018K^)7r, (15) 

^ I / 1-|(^t)iiP .... 

l(5'T)i2| = y ^ , (16) 

Arg (5t)i2 - - 0.153k + 0.585^2 - 0.087^^^ vr, (17) 

(5t)i3 = (5t)i2, (18) 

'^^^^^^'^ 1 + (1.780«: + 0.015«:T' ^^^^ 

Arg (5t)22 ^ (1 + 0.415/t + 0.508^2 + 0.089^^)77, (20) 

l(5T)23| = \/^^f|^-|(5T)22P, (21) 

Arg (5t)23 - + 0.205/t + 0.616^^ - 0.202/t3^ vr, (22) 

(5't)33 = (5't)22, and (23) 

^T* = ^T- (24) 



For K < 1, the above approximations for 15*1 and Arg {S) approximate the re- 
sults in Figs.El^c) and[2[d) respectively up to A\S\ < 0.01 and kig{S)] < 
0.01 accuracy. The exact equalities are due to the symmetry and unitarity of 
the S'-matrix, and the exchange of leads 2 and 3. 
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3 Comparison between reduced and unreduced systems 



In this section, tliree wave guides which are composites of the discussed L-bend 
and T-junction are analyzed. The scattering amphtudes from the original 2D 
structures and the reduced multiply-connected ID structures are compared. 
Cases of far apart and close discontinuities, different orientations of guides, 
and straight and smoothly curved guides are considered. 

The first example is a 2D square loop resonator with two leads as depicted 
in Fig. [3](a). The translational invariant sections have the same width W. For 
simplicity, the distances between the discontinuities are chosen to be the same, 
and are denoted by d as shown. In Fig.[3]^b), a reduced version for the structure 
in Fig. El^a) is shown. The magnitude of the transmission scattering amplitude 
15121 has been plotted versus k for the unreduced and reduced systems, in 
Figs. [3](c) and[3]^d), for d = 4W and d = O.IW respectively. The 5*12 for the 
2D structure is the scattering amplitude from the first transverse mode in lead 
2 to the first transverse mode in lead 1, and the dimensionless wave number k 
is now defined by k = k^^^W/n or k = kW/ir depending on the context. The 
procedure for a full- wave evaluation of |S'i2| for the 2D structure is standard 
[2Tp2] and it is not to be repeated here, but only the result is given. The result 
here is obtained with a FD TISE, and the number of sites across a width W 
is equal to 20. 

In the reduced system shown in Fig. [3t^b), a line labeled by rj and given a coor- 
dinate Xr^ has a wave function in the form given in Eq. ([2]). For the coordinates 
defined in Fig. [3](b), the wave functions on the lines are connected by 
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and 



1 + 
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This contains 10 equations with 10 unknowns, when Ai and A2 are given. The 
scattering amphtude S12 is obtained as 5*12 = -Bi at Ai = and A2 = 1. The 
magnitude 15*121 is plotted in Figs. [3](c) and[3]^d), for values of corresponding 
to the original 2D structure. In the calculation, Si^ and S't use the numerical 
values shown in Fig. [21 



In the case of d = 4W [see Fig. [3]J^c)], results from the original 2D and the 
reduced systems are nearly indistinguishable when seen in the size of the plot, 
while in the case of ci = O.IW [see Fig. [31(d)], the two results have a perceptible 
difference, especially when k becomes large. A criterion for the applicability 
of the one-mode reduction scheme is |/c*-^''(i| 3> 1, which means that the expo- 
nential tails of the evanescent waves for the higher transverse modes emanat- 
ing from the discontinuities, is much shorter than the distances between the 
discontinuities. In this regime, the evanescent waves from neighboring discon- 
tinuities do not overlap each other, and the discontinuities communicate with 
each other only via the first transverse modes in the translational invariant 
sections. Hence, the higher transverse modes in the sections are redundant, 
and a single-mode description is adequate. 

The length of an exponential tail of the second transverse mode is of the 
order of The values of are approximately equal to 0.18 

{k = 0), 0.23 {k = 1.0), 0.28 {k = 1.3), 0.37 {k = 1.5), 0.48 {k = 1.6), 
0.96 {k = 1.7), and 1.56 (k = 1.72), for the values of k given in the brackets 
[23j . It is seen that becomes large only when k approaches v^- For 

d = 4:W, \k^'^^d\ ~ 14 ^ 1 at K = 1.3, and that justifies the reduction scheme 
in the entire range of k in Fig. [31(c) [and also Fig. [D^c) later in this section]. 
For d = O.IW, \k^'^^d\ ~ 0.54 at k = 0, and \k^'^^d\ is certainly not "large" for 
Fig. [3]^d) [and also Fig. [D^d) later in this section] . In spite of this, the scheme 
might still perform up to certain precision, until it really starts to breakdown 
at ft ~ 1.2 as seen in Fig. M^d) [see also Fig. [11(d)]. However, its reliability in 
this regime of k is uncontrolled in general. 

Another 2D wave guide as shown in Fig. [D^a) is also analyzed. This struc- 
ture resembles the one in Fig. [3]^a), except that one of the leads is rotated 
by 90-degree. The width of the uniform sections are also denoted by W, and 
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the distances between the discontinuities are also denoted by d. A reduced 
structure for the guide is shown in Fig. \Mjo), and the magnitude of the trans- 
mission scattering amphtude |S'i2| is also plotted for both of the reduced and 
unreduced structures in Figs. |l](c) [for d = AW] and|l](d) [for d = O.IW]. 

The connecting equations for the reduced structure in Fig. m^b) are Eqs. (^S^, 
(ESD, (EZD, and 



1 + ^1 
ik 
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(29) 



The transmission scattering amplitude 5*12 is found as in the previous example. 
Likewise, it is seen in Fig. Hl^c) that the reduction scheme is guaranteed to work 
in the |A;*^^''(i| ^ 1 regime. Also, it is seen in Fig. ID^d) that the scheme might 
still work qualitatively or semi-quantitatively, when k is departed from this 
regime. 

Comparing the result in Fig. [3]^c) with that in Fig. ID^c), and the result in 
Fig. [3](d) with that in Fig. 111(d), it is seen that scattering amplitudes can 
depend significantly on the orientations of the branch guides at a discontinuity. 
This indicates that reduction schemes with symmetric branch lines such as 
those in Refs. [T]|2][51[T7] are not adequate for some cases. 

The third example is a 2D annulus structure with an inner and an outer 
radii of i? — W jl and R + W jl respectively, and two mutually perpendicular 
leads of width W radially connected to the annulus as shown in Fig. [S](a). 
Figure [5]^b) shows a reduced version of it. For the 2D annulus, we may follow 
a mode- matching full- wave treatment formulated by Xia and Li [2l]. In the 
annulus the wave function \l/ann. can be expanded by radial and angular modes, 
\E'ann.(^, 6') = Y^fl-M '^li^''")^^''^ ^ wherc a radial mode is given by (f)i{Kr) = 
aiJi{Kr) + biYi{Kr), and K = ^/2moE/h. The r and 9 are the radial and 
angular coordinates respectively; and the Ji and Yi are the Bessel functions of 
the first and second kinds respectively. At the inner radius, 4>i\r=R-w/2 = for 
any 9; At the outer radius, '^a.nn.\r=R+w/2 = when 9 is away from the leads, 
and ^ann. |r=i?+Ty/2 = wheu 9 is in the range of lead rj. In addition, the radial 
derivative d'ifi^rm./dr is equated with the longitudinal derivative when 
^ann. and meet at the outer arc of the annulus. The difference between 
the straight transverse cuts of the leads and the outer arcs of the annulus is 
neglected. The wave functions in the leads and the annulus are hence matched, 
and one can get a set of equations relating the coefficients of the modes in 
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the different regions. The transmission scattering amphtude 1 5*121 for the 2D 
annulus is plotted for R = ?>W [Fig. EKc)] and R = 0.8W [Fig. E^d)]. 

To apply a reduced calculation to the reduced system in Fig. Mjo), note that 
for smoothly curved guides with small curvatures, the back-scattering is small 
and the guides can be treated as refiectionless for most purposes. Guides 
with constant curvatures are translational invariant and indeed refiectionless, 
though the lengths and widths may not be rigorously defined. The junctions 
are treated as T-junctions. For the left junction. 
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ik 
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for the right junction. 
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(30) 
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(31) 



We let dr^ = R9r, — W, where 6*3 = 37r/2 and 64^ = Tr/2, and the subtraction is to 
approximately exclude the regions in the junctions [see region D in Fig. [D^b)]. 
The radius R = 3W gives ds ~ 13.1W and d^ ~ 3.7W, and R = 0.8W gives 
ds ~ 2.77W and d^ ^ 0.26W. Results for |S'i2| are plotted in Figs. [5](c) and 
Eld). 

If one analyzes the reduced system using the Griffith scheme [Tp|3f6] for a 
junction with three branches, the left junction has 
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and the right junction has 



dxA 



(32) 
(33) 



X4=d4 



2\x2=0 



3 |x3=d3 



412:4=0 



and 



(34) 



10 



+ 



di/j4 



0. 



(35) 



dx2 



X2=0 



X3=d3 



dX4 



X4=0 



Note that the equations at a junction are symmetric with respect to an inter- 
change of any two branches. Using the same ^3 and d^, results for |S'i2| are 
also plotted in Figs. [5](c) and[5]^d). Comparing the results from our scheme 
and the Griffith scheme with the result from the 2D calculation, it is seen that 
the Griffith result is qualitatively different from the 2D result in general, while 
our scheme captures the essential features in the 2D result, especially in the 
case of longer guide sections [Fig. [5](c)]. 

The lengths of the guide sections in Fig. [S] are about the sizes of those in 
Figs. [3] and IH but the disagreement between the 2D and reduced calculation 
results are seen to be more severe in Fig. [5l especially in the case of short 
guide sections [Fig. [5]^d)]. This is mainly due to vaguer notions of the lengths 
and widths of the guide sections between the discontinuities, and a stronger 
distortion of the shapes of the junctions from a "T" [Fig. [Il^b)], when the 
guide sections are curved and short, in addition to a violation of the criterion 
^ 1. Comparing Figs. [5](c) and[5]|^d), it is also seen that our scheme 
performs better when the discontinuities are more apart. 



4 Concluding Remarks 



This study has shown that the higher transverse modes in the uniform-cross- 
sectional sections in wave guides can give only a minor effect on the low-energy 
properties of the guides. A scheme to remove these modes and a criterion for 
the energy range ^ 1) in which the scheme applies have been proposed. 

In the scheme, a reduced system and its corresponding system in a realistic 
space have a sound relationship, and therefore the scheme may find more 
practical use in the analyses of realistic wave guides than previously reported 
schemes p!|17fl8j . 

The electronic spin degrees of freedom may also be included into the scheme by 
expanding the current one-mode S'-matrices to two-mode S'-matrices, for the 
up and down spin channels, and in the same way as in Sec. [3l the two-mode S'- 
matrices are calibrated by full-wave calculations for systems in realistic spaces. 
Likewise, the precision of the scheme at short guide sections can also be refined 
by including more transverse modes and using multi-mode S'-matrices. 



11 



References 

[I] H. Kuhn, Helv. Chim. Acta 32, 2247 (1949). 

[2] J. Stanley Griffith, Trans. Faraday Soc. 49, 345 (1953); ibid., 49, 650 (1953). 

[3] K. Ruedenberg and C. W. Scherr, J. Chem. Phys. 21, 1565 (1953). 

[4] P. Exner and P. Seba, Rep. Math. Phys. 28, 7 (1989). 

[5] T. Kottos and U. Smilansky, Ann. of Phys. 274, 76 (1999). 

[6] J.-B. Xia, Phys. Rev. B 45, 3593 (1992). 

[7] J. M. Mao, Y. Huang, and J. M. Zhou, J. Appl. Phys. 73, 1853 (1993). 
[8] P. Singha Deo and A. M. Jayannavar, Phys. Rev. B 50,11629 (1994). 
[9] M. V. Moskalets, Low Temp. Phys. 23, 824 (1997). 
[10] C.-M. Ryu and S. Y. Cho, Phys. Rev. B 58, 3572 (1998). 

[II] S. Bandopadhyay, P. Singha Deo, and A. M. Jayannavar, Phys. Rev. B 70, 
75315 (2004). 

[12] D. Bercioux, M. Governale, V. Cataudella, and V. M. Ramagha, Phys. Rev. 
Lett. 93, 56802 (2004). 

[13] P. Foldi, B. Molnar, M. G. Benedict, and F. M. Peeters, Phys. Rev. B 71, 33309 

(2005) . 

[14] U. Aeberhand, K. Wakabayashi, and M. Sigrist, Phys. Rev. B 72, 75328 (2005). 

[15] X. F. Wang and P. Vasilopoulos, Phys. Rev. B 72, 165336 (2005). 

[16] O. Kalman, P. Foldi, M. G. Benedict, and F. M. Peeters, larXiv:0806.2734] 
(unpublished). 

[17] K.-K. Voo, S.-C. Chen, C.-S. Tang, and C.-S. Chu, Phys. Rev. B 73, 35307 

(2006) . 

[18] B. Shapiro, Phys. Rev. Lett. 50, 747 (1983). 

[19] M. Buttiker, Y. Imry, and M. Y. Azbel, Phys. Rev. A 30, 1982 (1984). 

[20] Y. Gefen, Y. Imry, and M. Y. Azbel, Phys. Rev. Lett. 52, 129 (1984). 

[21] S. Datta, Electronic Transport in Mesoscopic Systems, 1st Ed. (Cambridge 
University Press, Cambridge, 1995). 

[22] D. K. Ferry and S. M. Goodnick, Transport in Nano structures, 1st Ed. 
(Cambridge University Press, Cambridge, 1997). 

[23] /c(2) and k are related by k^'^^W = tt^J - 3. 

[24] J.-B. Xia and S.-S. Li, Phys. Rev. B 66, 35311 (2002). 



12 



(a) 



(b) 



(c) 




x=x = 

1 2 



(d) 2 3 



X=X=X=Q 



1 



Fig. 1. (Color online) The two discontinuities in 2D wave guides considered in this 
paper, (a) the L-bend and (b) the T-junction. The widths of the guides are denoted 
by W . The branch guides are labeled by 1, 2, and 3, in which coordinates are defined 
as shown. The translational noninvariant regions are labeled by D. In the reduced 
systems, (c) a L-bend is represented by a "L," and (d) a T-junction is represented 
by a "T." The arrows on the branch lines indicate the positive directions of the 
coordinates defined on the lines, and the origins of the coordinates are defined to 
be at the branching nodes. 
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Fig. 2. (Color online) The (a) magnitudes of the elements of S'l, (b) arguments 
of the elements of 5l, (c) magnitudes of the elements of S't, and (d) arguments of 
the elements of S't are plotted versus a dimensionless longitudinal wave number k, 
defined by k = k'^^^Wj-K. 
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Fig. 3. (Color online) (a) A considered 2D composite structure. The uniform-cross- 
sectional sections have the same width W . The distances between the discontinuities 
are chosen to be the same and are denoted by d. (b) A reduced version of the above 
structure. A coordinate Xr^ is defined on line (r/ = 1 — 6) with a positive direction 
indicated by an arrow on the line. The coordinates of the discontinuities on the 
lines are shown in the boxes, (c) For d = AW , the magnitude of the transmission 
scattering amplitude |5'i2| is plotted versus k, for the original 2D system (solid line) 
and the reduced system (dotted line). The k here is defined by k = k^^^W/ir or 
K = kW/n. Seen in the size of the present graph, the two curves are almost indis- 
tinguishable, (d) For d = 0.1 VF, results for |5i2| are also shown. The curves deviate 
from each other, especially when k gets^^rger. 
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Fig. 5. (Color online) (a) A considered 2D annulus with two mutually perpendicular 
and radially connected leads. The leads have the same width W. The inner radius 
R — W/2 and outer radius R + W/2 of the annulus are defined as shown, (b) A 
reduced version of the annulus. A coordinate Xjj is defined on line = 1 — 4) 
with a positive direction indicated by an arrow on the line. The coordinates of the 
junctions on the lines are shown in the boxes, (c) For R = 3W, the magnitude 
of the transmission scattering amplitude 15*121 is plotted versus k, for the original 
2D annulus (solid line) and the reduced system (dotted line for our scheme, and 
dash-dotted line for the Griffith scheme). The k here is defined by k = k^^^W/ir 
(note that k^^^ is defined in the leads) or k = kW/ir. (d) For R = 0.8VF, results for 
|5i2| are also shown. 



17 



